Various spin-polarization states beyond the maximum-density droplet: 

a quantum Monte Carlo study 
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Using variational quantum Monte Carlo method, the effect of Landau-level mixing on the lowest- 
energy-state diagram of small quantum dots is studied in the magnetic field range where the density 
of magnetic flux quanta just exceeds the density of electrons. An accurate analytical many-body 
wave function is constructed for various angular momentum and spin states in the lowest Landau 
level, and Landau-level mixing is then introduced using a Jastrow factor. The effect of higher 
Landau levels is shown to be significant; the transition lines are shifted considerably towards higher 
values of magnetic field and certain lowest-energy states vanish altogether. 
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Quantum dots (QD) are systems containing a num- 
ber of charge carriers in a nanoscale volume. A two- 
dimensional semiconductor QD |jj can be constructed 
e.g. by superimposing a confining potential on two- 
dimensional electron gas, which can be fabricated into 
the inversion layer between two different semiconductor 
materials. In the inversion layer, the confinement per- 
pendicular to the interface can be made so strong that 
only the lowest-energy eigenstate in that direction has a 
nonnegligible probability of being occupied. This makes 
the system genuinely two-dimensional, and leads to some 
physical effects that are absent in higher dimensions: in 
thcrmodynamical limit the integer and fractional quan- 
tum Hall effects Q are perhaps the most famous exam- 
ples. 

Experimentally, QDs are observed to have a series of 
different ground states as the magnetic field is increased 
H . While many properties of the states can be measured, 
understanding their nature is still one of major theoret- 
ical goals in the field. In a large range of moderately 
high magnetic field the very stable and fully spin polar- 
ized maximum-density-droplet state (MDD) remains the 
ground state. This many-particle state has an extremely 
accurate analytical description [Q, and is here used as 
the starting point to study phase diagrams of possible 
lowest-energy states in stronger fields. 

The variational quantum Monte Carlo (VMC) || 
method is used to study the QD system. In this method 
a r trial many-body wave function ^({i^}, with de- 

sired properties and with free variational parameters on 
is first constructed, and then the parameters are opti- 
mized to converge towards the exact wave function ho- 
using the optimized wave function, the expectation value 



of an observable A can be evaluated as the average of the 
corresponding local quantity ^^ X A^, e.g., for the Hamil- 
tonian operator Ti: 
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(1) 



where the A^-particle-coordinate configurations are 
distributed as \§>\ 2 and generated using the Metropolis 
algorithm. 

The variational principle guarantees that the total en- 
ergy given by the VMC method, using any trial wave 
function, is always an upper bound for the true total 
energy of the quantum state in question. The variance 
of the local energy \& _1 7^5' diminishes as the trial wave 
function approaches an eigenstate of the Hamiltonian, 
and as a result it can be used not only as a measure of 
the statistical error in E^, but also as a measure of the 
difference between calculated and true energies Eq,—E^ . 

The variational parameters in the trial wave function 
are optimized by minimizing the total energy. The mini- 
mization process itself was done using the stochastic gra- 
dient method [[| with analytical expressions for deriva- 
tives . The method has proven to be fast in finding a 
minimum of a function whose values are not exact. 

The QD in this study is an A^-electron system on a 
two-dimensional plane, in rotationally symmetric and 



parabolic potential V(r) 



and in perpendicu- 



lar magnetic flux density B = Bu z . An effective mass m* 
is used to describe the effects of the underlying crystal 
structure, and o>o determines the strength of the in-plane 
confinement. In the symmetric gauge, A = x r, the 
effect of the vector potential on the Hamiltonian is to en- 
hance the confinement strength, cj 2 — > Wq + (cj c /2) 2 = cj 2 , 
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where lo c — eB/m* is the effective cyclotron frequency, 
and to introduce angular-momentum and spin-dependent 
terms. Switching to the effective harmonic oscillator 
units (to* = h = u) = e = 1), the total Hamiltonian 
operator can be written as 

N r 

h = 3 E t- v ? + r " + *s*,< + 7x0] + E - ■ 

Z— 1 1<J J 

(2) 

Here l z ,i and s 2) j are the z-components of angular- 
momentum and spin operators for the ith particle, 7* is 
defined as 7* = g*m* /mo, where g* is the effective Lande 
g-factor, and C is a dimensionless interaction strength, 
C 2 = Ha,/tko ■ m*/(e 2 m ). In these units lu c = B < 2. 

Single-particle states of the noninteracting part of the 
Hamiltonian (|^) can be solved analytically for arbitrary 
magnetic flux density B §. As B is increased the states 
separate to Landau levels with level spacing asymptoti- 
cally proportional to B, and intra-level energy spacing of 
states asymptotically proportional to B~ x . In the limit 
of infinite magnetic flux density only the lowest Landau 
level (LLL) remains occupied. The single-particle wave 
functions (p m , that the LLL is composed of, have a par- 
ticularly simple form: 

^ m (z) = zl" l lexp(-i|z| 2 ), (3) 

with energies (in units of Tilo) 

E m , = \ + {\-\u3 c )m. (4) 

In Eq. (^) , z is a complex coordinate in the plane of the 
electrons: z = x + iy, and m = 0, 1, . . . is l z /%. 

It is a common approximation to truncate the full 
Hilbert space to the lowest Landau level only, and ignore 
the effects of higher Landau levels altogether. For exam- 
ple, exact diagonalization calculations are done within 
the LLL approximation for N > 4. The approximation 
improves with increasing B, since higher Landau levels 
move farther off in energy, but as will be shown, Landau- 
level mixing (LLM) can qualitatively alter the behavior 
of the system even in such regions of B, where the LLL 
approximation is frequently used. 

In this study, the trial wave functions '5 are of the 
Jastrow-Slater form. The construction of if? begins 
by creating an LLL-many-particle wave function 
which determines, e.g., the angular momentum and the 
spin of the state. The LLM is then introduced using 
two-body correlation functions of Jastrow type [0]. To 
construct the LLL part of the wave function, the unnor- 
malized maximum-density droplet state $mdd is used as 
a starting point: 

N 

*mdd = \\{z 3 - Zi) exp(-±^ \ Zl \ 2 ) . (5) 

i<j i—1 



In the thermodynamical limit (luq — > 0) this configura- 
tion corresponds to the very stable quantum-Hall state at 
filling factor v = N/N& = 1 @. For QDs in the LLL ap- 
proximation, $mdd is the lowest-energy state in a large 
region around v = £mdd / L = 1 . The state is composed 
of consecutive single-particle states of the form of Eq. 
(||) with to = 0, . . . , N — 1, and it has the total angular 
momentum £mdd = ^N(N — 1). 

In this study only post-MDD states are considered, 
that is, states that have nonnegative additional angular 
momentum AM = L — Lmdd with respect to the MDD 
state. Since the parabolic form of the external confine- 
ment allows one to perform the elimination of the center- 
of-mass motion, the spin-up and spin-down degrees of 
freedom are not independent. Therefore it is assumed 
that it is favorable for an electron to put zeros of the 
wave function on the positions of other spin type elec- 
trons also. This is equivalent to the requirement that 
any post-MDD state should contain the ^mdd as a fac- 
tor: 

*LLL = *MDD PAM , (6) 

where Pam is a polynomial of proper symmetry and de- 
gree AM. Under this assumption, constructing the LLL 
trial wave function is a matter of choosing the suitable 
polynomial Pam- It must have the correct symmetry, 
and its degree AM gives the additional angular momen- 
tum with respect to the MDD state. 

The polynomial Pam is constructed as follows: One 
starts from the product 

AM 

Pam = II * ' ( ? ) 

»=i 

where the coordinate transformation £j = Zi — Zj 
has been applied to remove the center-of-mass motion. 
This fixes the additional angular momentum AM of the 
trial state. Correct symmetry is now built in using the 
Young symmetrization operator y x : 

*LLL = y x (*MDD Pam) = *MDD Y x P A M , (8) 

where y' x P& M is Pam above. The parameter x is the 
number of inverted spins with respect to a fully spin- 
polarized state, and it determines the shape of the Young 
tableau corresponding to y x . Because an antisymmetric 
factor $mdd was taken out of y x , the directions of sym- 
metrization and antisymmetrization are interchanged in 
y' x . By fixing the spin and angular momentum through 
parameters AM and Xt the LLL part of the trial wave 
function is determined completely, i.e., it has no varia- 
tional parameters. This family of trial wave functions 
covers most spin configurations in the region AM = 
0, . . . N. For example, for seven electrons there are only 
three states (out of 23 nondegenerate states) that are un- 
available: S = 1/2 at L = 27, and S = 3/2 and S =1/2 
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FIG. 1: Phase diagrams for six electrons in the lowest Landau level approximation (left panel) and with Landau-level mixing 
(right panel). The labeling of the states is: Nf, AM, where Nf is the number of spin-up electrons and AM is the additional 
angular momentum (Imdd = 15 for six electrons). The vertical axis is the strength of the Zeeman coupling per spin, AE Z = 
||/iBff*[, the value of which in GaAs (12.7/ieV/T) is marked by dashed lines in the figures. Other parameters were: m* /mo = 
0.067, Tiuia = 5meV and e r = 12.4. The relative interaction strength C varies from 1.23 {B = 7T) to 1.01 (B = 12 T). 



at L = 28. These inaccessible configurations are not ex- 
pected to be lowest-energy states in any case. For states 
with 2S — \N — 2 AM |, the LLL construction above re- 
sults in wave functions that are very similar to skyrmion 
states p"o| . 

After the lengthy construction of the LLL part, the 
LLM part is simply multiplied into the trial wave function 
as a two-body correlation function of Jastrow type 



J\\ r ij}i<j) 



JJexp 

i<3 



1 + aan 



The total wave function is then 



* = #LLL J ■ 



(9) 



(10) 



This form of the correlation factor leaves the spin and 
angular momentum properties of \& intact. The Jastrow 
factor has two adjustable parameters: ay = for par- 
allel spins and ay = a-fj. for antiparallel spins. For each 
trial state, these parameters are functions of B, and have 
to be optimized separately for each magnetic-field value. 
Constants /?y are determined by cusp conditions [ pd| . 
The form of Eq. (|l0|) has proven very efficient in captur- 
ing the Landau-level mixing; in small systems as much 



98% of the correlation energy can be recovered 1 12 

Results for two systems, N — 6 and N = 7, are pre- 
sented in this article. In order to test the accuracy of the 
LLL construction an exact diagonalization calculation 
p3[ was performed for seven electrons at B = 5 T. This 
showed that for all possible LLL lowest-energy states the 
error in total energy is very small, at most 0.065%. In the 
case of the LLM wave function, no such diagonalizations 
are currently available. 



Fig. [j] shows results for a six-electron QD, and corre- 
sponding data for seven electrons are shown in Fig. |^. 
Both systems use GaAs parameters: m IC \ = 0.067, 
g * = -0.44, and e r = 13.0 for N = 7 as in but 
12.4 for N = 6 as in |4j. 

The figures show magnetic flux density vs. Zeeman 
coupling strength AE Z = ||/ib.9*|- Generally, with in- 
creasing Zeeman coupling strength the degree of spin po- 
larization increases, and above some critical value AE* 
the system remains spin-polarized at all values of B in the 
post-MDD region. In the LLL approximation, this hap- 
pens at AE* Z = 26.5 ^eV/T for N = 6, and at 34.5 /xeV/T 
for N = 7. Lowering the Zeeman coupling strength al- 
lows other total-spin values, and with AE Z — the sys- 
tem goes through all possible spin configurations. The 
lowest-energy states are the skyrmionic states as dis- 
cussed above. 

The effect of the Landau-level mixing can be seen 
on the right-hand panels. Quantitatively, the transition 
lines are shifted towards higher fields and lower Zeeman 
couplings. In the case N = 6, half of the lowest-energy 
states even vanish completely when the LLM is taken 
into account. 

The shift towards higher B means that in the LLM 
case the states are more stable against radial expansion 
as B is increased. (AM increases as one moves to the 
right in the figures.) This should be compared with Rcf. 
|l5| , where similar results have been obtained within the 
LLL approximation, by assuming the QD to have finite 
thickness. Both effects are results of the expansion of 
the basis set of the Hilbert space, but in our case the sys- 
tem still remains two-dimensional. The vanishing of some 
lowest-energy states can be explained by noting that the 
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FIG. 2: Phase diagrams for seven electrons. In this case Lmdd = 21, and the value e r = 13.0 was used to compare with an 
exact diagonalization calculation. The relative interaction strength C varies from 1.38 (B = 3T) to 1.17 (B = 7.0 T). 



two-body correlation factor has strongest effect on the 
states that are most compact (i.e., have the lowest an- 
gular momenta). For the MDD state, the gain in energy 
is largest, and it can completely block some neighbor- 
ing states whose energy is only slightly lower in the LLL 
approximation. 

The value of AE Z (corresponding to g* = —0.44) for 
GaAs is 12.7/xeV/T. This is a maximum value (within 
approximations and assumptions used) that can only be 
lowered, using tilted field experiments. Since the value is 
much lower than AE* , it should be possible to encounter 
partially spin-polarized post-MDD states in experiments, 
even without tilting the magnetic field. 

In summary, simple trial wave functions for partially 
and fully spin-polarized QD systems have been con- 
structed. The wave functions were shown to produce 
excellent total energies, and provide an accurate estimate 
of Landau-level mixing. 

We demonstrate that the LLM, which is almost al- 
ways neglected in previous studies, is able to suppress 
the existence of certain lowest-energy states with small 
amount of additional angular momentum (compared to 
N). Furthermore, a strong shift of transition points to- 
wards higher magnetic fields due to the LLM is observed. 
Another important result is, that despite the LLM, par- 
tially spin-polarized states can exist in the post-MDD re- 
gion with realistic Zeeman coupling strengths. The most 
visible of the these states is likely to have a single spin 
flipped and angular momentum equal to Lmdd + N — 1. 
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(2000-2005). 
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